Computation of spectroscopic factors with the coupled-cluster method 
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We present a calculation of spectroscopic factors within coupled-cluster theory. Our derivation 
of algebraic equations for the one-body overlap functions are based on coupled-cluster equation-of- 
motion solutions for the ground and excited states of the doubly magic nucleus with mass number 
A and the odd-mass neighbor with mass A — 1. As a proof-of-principle calculation, we consider 
16 and the odd neighbors 15 and 15 N, and compute the spectroscopic factor for nucleon removal 
from 16 0. We employ a renormalized low-momentum interaction of the V[ w-fc type derived from a 
chiral interaction at next-to-next-to-next-to- leading order. We study the sensitivity of our results 
by variation of the momentum cutoff, and then discuss the treatment of the center of mass. 
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I. INTRODUCTION 



In the past two decades, ab initio nuclear structure 
calculations have led to the development and test of 
high-precision models with predictive power [H-Ilj]. Re- 
cently, the application of effective field theory (EFT) 0- 
0] and renormalization group techniques [H, |9| resulted in 
a mo del- independent approach to the nuclear interaction. 
These approaches have significantly deepened our under- 
standing of nuclear forces and have also provided us with 
new technical means to simplify the solution of the nu- 
clear many-body problem. The interactions from chiral 
EFT have been probed in light nuclei flfj| - fl3j and selected 
medium- mass nuclei with different techniques [l4l - fl6j ]. 
The focus of ab initio calculations is not only on observ- 
ables such as binding energies, radii, and low- lying exci- 
tation spectra, but also on transition rates and more de- 
tailed spectroscopic information. Very recently, ab initio 
theory be gan to bridge the gap from nuclear structure to 
reactions [17H19I ]. The inclusion of continuum effects, for 
instance, is necessary for the description of weakly bound 
and unbound nuclei. Direct reactions such as stripping 
and pickup of a single nucleon are rather well understood 
within phenomenological approaches (see, e.g. Ref. (20|). 
but constitute a current frontier for ab initio theory. 

The interpretation of direct reactions within a given 
model or Hamiltonian is based on spectroscopic fac- 
tors [2l|, [22j ■ The spectroscopic factor depends on wave 
function overlaps (see Eq. (TJ5| below for a definition) and 
provides useful information that relates nuclear struc- 
ture within a given model (i.e. within a given Hamilto- 
nian) to stripping and transfer reactions [2l|. The spec- 
troscopic factor is not an observable as it depends on 
the employed Hamiltonian or model. In nuclear physics, 
the high-momentum parts of the interaction are uncon- 
strained and modeled in different ways. Thus, the short- 
ranged part of the wave function is mo del- dependent, and 
so is an overlap between wave functions. Therefore, the 
spectroscopic factor is merely a theoretical quantity and 
cannot be measured [23|,[24|]. However, the spectroscopic 



factor "provides a useful basis for the comparison of ex- 
periment and current nuclear models" [2l| . Its purpose 
thus lies in understanding a direct reaction within a cer- 
tain model or Hamiltoniam and this interpretation might 
be useful and interesting (251 ] . 

In this paper, we develop the technical tools to com- 
pute spectroscopic factors within the coupled-cluster 
method [2g-[3fJ (see Ref. [31( for a recent review of this 
method), and perform a proof-of-principle calculation for 
16 O. The computation of the spectroscopic factor within 
coupled-cluster theory is not trivial (i) since the method 
does not readily yield a many-body wave function, and 
(ii) due to details related to the translation invariance of 
the coupled-cluster wave function. This paper is struc- 
tured as follows. Section [H] is dedicated to a summary 
of the employed coupled-cluster method. The theoreti- 
cal computation of spectroscopic factors within coupled- 
cluster theory is presented in Section Hill We present our 
results and a discussion of the center-of-mass treatment 
in Section ITVl Section IVl contains our conclusions and an 
outlook. 



II. EQUATION-OF-MOTION AND 
COUPLED-CLUSTER THEORY FOR NUCLEI 

In this section we introduce the Hamiltonian and 
coupled-cluster theory (26l - [3fj| for closed-shell and open- 
shell nuclei. Although our implementation of coupled- 
cluster theory has been presented elsewhere [HI, l3lU34j , 
we give a brief overview of the method since some details 
are needed for the calculation of spectroscopic factors. 

We consider the intrinsic nuclear A-body Hamiltonian 
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Here T is the kinetic energy, T cm is the kinetic energy 
of the center-of-mass coordinate, and V is the two-body 
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nucleon-nucleon interaction. In this paper we use low- 
momentum interactions Vjow-fc foL f35jj with sharp cutoffs 
A = 1.6,1.8,2.0,2.2 fm _1 , respectively. For simplicity, 
we neglect any contributions of three-nucleon forces as 
we focus on a proof-of-principle calculation. 

In coupled-cluster theory, one writes the ground-state 
many-body wave function as 

|V>o) = e T |0o) ■ (2) 

Here, \<po) is a product state. The cluster operator T in- 
troduces correlations as a linear combination of particle- 
hole excitations 

T = T 1 +T 2 + ... + T A . (3) 
Here, the n-particle-n-hole excitation operator is 
/ 1 \ 2 " 

T « = b ) II E ■ ■ ■ -s- (4) 

f = 1 Civ ,iu 

We employ the standard convention that indices ijk . . . 
refer to orbits below Fermi level (holes) and abc . . . above 
Fermi level (particles). Approximations in coupled- 
cluster theory are introduced by truncating the cluster 
operator T at a certain particle-hole excitation level. In 
this work we truncate T at the two-particle-two-hole 
excitation level, i.e. T m T\ + T2, which gives the 
coupled-cluster method with singles and doubles exci- 
tations (CCSD). This is the most commonly used ap- 
proximation, as it provides a good compromise between 
computational cost on the one hand and accuracy on the 
other. 

Within the CCSD approximation, the computational 
cost is given by r? n\, where n and n u denote the number 
of occupied and unoccupied orbitals, respectively. The 
correlated ground state solution is given by the ampli- 
tudes t\ and tfj that solve the nonlinear equations 

(#|2F|to) = , (5) 

«|#|4>o> = . (6) 

Here, the bra states are particle-hole excitations of 
the reference Slater determinant, and H denotes the 
similarity-transformed Hamiltonian, 

H = e- T He T =(He T ) c . (7) 

The subscript c indicates that only fully connected dia- 
grams give non-zero contributions. Once T is determined 
from the solution of the coupled-cluster equations (J5j , the 
correlated ground state (g.s.) energy is given by 

E = ((fo\H\<h) ■ (8) 

The CCSD approach is known to work particularly well 
for the ground state of nuclei with closed (sub-) shells, as 
a Slater determinant provides a reasonable first approx- 
imation. In this work we use the Equation-of-Motion 



(EOM) [H IH-dil method to solve for the ground and 
excited states of the closed-shell nucleus A and its odd 
neighbors with mass number A — 1. 

In EOM, the ground and excited states of a nucleus 
with mass number B are obtained by acting with an ex- 
citation operator £2^ on the ground state wave function 
of a nucleus with mass number A, i.e. ip^ — fl^tpQ. 
Here \i denotes quantum numbers such as spin, par- 
ity, and isospin projection. Within the EOM approach, 
the ground state wave function tp A denotes the coupled- 
cluster wave function e T (f>Q. In this work we choose either 
B = A, in which case we solve the excited states of closed- 
shell nucleus A, or B = A — 1, in which case we solve the 
ground and excited states of the A — 1 neighboring nu- 
cleus. To solve for the excited states of the closed-shell 
nucleus A, we define £2 M by the excitation operators, 

R A = r + J2 rfalat + ± £ r^a^-a, , (9) 

i a ijab 

l a = 1 + W<>a + ~ E &«W a »°» ■ ( 10 ) 

ia ijab 

Here, we suppressed the index fi, but it is understood that 
the operators R A and L A excite and de-excite states with 
quantum numbers fx, respectively. For the ground and 
excited states of the nucleus with mass number .4—1, 
we define £2^ by the particle removal operators 

r*- 1 = J2 r * a * + 2 E r ?A a o a i > ( n ) 

i ija 

L A ~i = Y,M + \Y.?i44a a . (12) 

i ija 

Again, we supressed the index \x labeling the quantum 
numbers. The operators i? M = R A (R^ = R A ~~ 1 ) com- 
mute with the cluster operator T, and the unknowns 
rf , rfj , r°- ) solve the EOM equation 

pT.flJ \4> ) =u>M(h) , (13) 

which defines an eigenvalue problem for the excitation 
operator with eigenvalue ui^ = — Eq. It is clear 
from the definitions ((H) and ([7]) that H is non-Hermitian; 
dual space solutions need to be calculated explicitly. We 
obtain the de-excitation operators L M = L A ,L A ~ X by 
solving the left eigenvalue problem 

(<l> \L fi H = (0 O |V^ • (14) 

The right and left eigenvectors form a bi-orthogonal set 
and are normalized in the following way 

(^olMVl^oH <W • (15) 

The EOM solution for the ground state of system A 
is identical to the CC solution, so that R A = r = 1. 
Ref [3l| provides a detailed description of EOM. 
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III. OVERLAP FUNCTIONS AND 
SPECTROSCOPIC FACTORS FROM 
COUPLED-CLUSTER THEORY 

The one-particle overlap function between two wave 
functions "5 a-i and ^a of nuclei with mass number A— 1 
and A, respectively, is defined as [22| 



Oi-i(2) 



(16) 



Here £ represent the 3 (A — 2) translationally invariant 
position coordinates and the A — 1 spin coordinates of 
A — 1 particles present in both \& a-i and ^a, while 2? 
labels the position and spin of the additional particle 
in the nucleus with mass number A with respect to the 
center of mass of the nucleus with mass number A — 1. 
The isospin coordinate has been suppressed. 

In our coupled-cluster approach, however, we do not 
employ coordinates with respect to the center of mass, 
as this would limit us to light systems [40|. Thus, the 
overlap can be associated with a specific nucleon repre- 
sented by a second quantization operator, 



Oi^W^iA-llaffilA) 



(17) 



Here, \A) and \A — 1) denote eigenstates in the nucleus 
with mass A and A — 1, respectively. Typically, \A) is 
the ground state, and \A — 1) is the ground state or an 
excited state. Our formalism will be kept general and is 
not limited to these cases. The radial overlap function 
A _i(lj; r) is derived by expanding x in terms of partial 
waves, 



a(x) = (-r™ 



E 

Ijm 



i(r)Yi jm (x) 



(18) 



Here, Yij m (x) is the spin-orbital spherical harmonic 

Y ljm (x) = \Yi{r)®Xi/2{a)] jm , 

Yiif) is the spherical harmonic of rank I and Xi/2( cr ) is 
a fermionic spin function. The orbital angular momen- 
tum quantum number is denoted by /, while j and m 
denote the rank and projection, respectively, of Yij m (x) 
as a spherical tensor. The hat denotes unit vectors, i.e. 
x = x/\x\. We have also introduced the spherical anni- 
hilation operator aij m (r) = (—) J+m aij- m {r). The radial 
overlap is now given by the reduced matrix element and 
the overlap becomes 



Oi_ x {x) 



E 

i 



{-) j - m (JaMaJ - m I Ja~iM A -i) 



xOi_.il j ]r )Y l3m (x) . (19) 

Here, (• | •) denotes a Clebsch-Gordan coefficient. The 
overlap is now expressed by a radial function associated 
with each tensorial component Yij m (x) 

Oi^ilj-r) = (A-l\\~ aij (r)\\A) (20) 
, m (A-lM A ^\a ljm (r)\AM A ) 



(-y 



This equation also defines the reduced matrix elements 
we employ. 

The norm of the radial overlap function is the spectro- 
scopic factor 

SU(lj) = J drr 2 \0U(ly,r)\ 2 . (21) 

The overlap functions can be expressed in an energy basis 
by inserting the expansions 

4™M = Y, a UmMr) , (22) 

n 

aij m (r) = ^a n i ]m (j)* nl:j (r) , (23) 

n 

where n is the nodal quantum number and 4> n ij (r) is the 
radial single-particle wave function associated with the 
orbits nljm. While a\- m (r) represents the creation of a 

particle at radial distance r, a n ij m represents the action 
of populating a single-particle orbit. 

Assuming orthogonality of the single-particle wave 
functions, the spectroscopic factor is written as, 



Si-im - ^2\{A-l\\~a nlj \\A) 



E 



I (A - l\a n ij m \A) | 2 
{J A M A j-m I Ja-iM^-i) 2 



(24) 



Here we used the Wigner-Eckhart theorem for the re- 
duced matrix elements. Eq. (|24|) is our starting point 
since we work in an uncoupled (m— scheme) basis. Using 
the EOM-CCSD solutions for the right and left eigen- 
value problems for the A and the A—l systems, and em- 
ploying the ground state solutions for system A, Eq. (|2"4")l 
takes the form, 



E 



(J A M A j-m I Ja-iMa-i) 2 



(25) 



This gives the equation for the spectroscopic factors as 
defined within coupled-cluster theory. We note that this 
equation is unambiguously and uniquely defined in terms 
of the left and right eigenstates of the nuclei with mass 
numbers A and A—l. This is clear since the spectroscopic 
factor is given by the absolute value squared of the one- 
body overlap matrix element, so any ambiguity related 
to the normalization condition (|15[) is removed. 

In Eq. (|2*5|) we have introduced the similarity- 
transformed creation and annihilation operators, 



— t t t 
e a^e 

-T T 

e a v e 



(26) 
(27) 



(J A M A j-m I J A -iM A _i) 



Using the Baker-Campbell-Hausdorff commutator ex- 
pansion, we can derive algebraic expressions for a\, and 
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A/ 


r 2 


V V 


L A 


♦ A ♦ A A 


A-l 


^+ \A 


R A 


♦ \L + V V 


A-l 


\ ♦ \ V 



TABLE I: The left panel lists algebraic symbols for the 
particle-hole excitation operators and particle removal oper- 
ators. Corresponding diagrammatic representations are dis- 
played on the right. Arrows pointing up (down) represent par- 
ticle (hole) orbits with implicit summation indices a,b,c, . . . 
(i.j,k,...). 

in terms of the "bare" creation and annihilation operators 
a p and a p and the particle-hole excitations amplitudes tf 



4 = ^-EW-^E^fe . ( 28 ) 

b jbc 

= a v+Yl iV i a% + \ ■ ( 29 ) 

i ijc 

These equations can also be given in a diagrammatic 
form, which provides a convenient bookkeeping system 
for the available Wick contractions in the expressions. 

The coupled-cluster diagrams are similar to Goldstonc 
diagrams. An algebraic Wick contraction corresponds to 
the diagrammatic connection of two directed lines, but 
the interpretation rules are slightly different. We refer 
the reader to Refs. jU, for a complete introduction 
to the diagrammatic approach. Here we will only present 
the few concepts necessary to introduce the novel exten- 
sions of the formalism used in the context of spectroscopic 
factors. 

Diagrammatic representations of the excitation and 
particle-removal operators T, L M , and are displayed in 
Table HI Lines with arrows pointing up (down) represent 
particle (hole) orbits. These lines have implicit indices 
a, b, c, . . . (i.j, k, . . .) that are summed over. We suppress 
both the summation symbol and the dummy indices for 
a cleaner notation. 

We have to deal with diagrams that represent opera- 
tors with an index that is not being summed over. Such 
a creation (annihilation) operator is represented by a di- 
rected line pointing out from (in to) a small circular ver- 
tex. The corresponding diagrams are displayed in the 
upper half of Table HI 

Eqs. ([25)) and can be reproduced diagrammatically 
as displayed in the lower half of Table |TTJ The possible 
Wick contractions between the creation and annihilation 
operators a) p and a p , and the cluster operators Ti and 
Ti depend on whether the index p denotes an orbital 
above or below the Fermi surface. The small circular 
vertices distinguish the index fixed by the operator and 
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ai, a a 


\ / 
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\ + V + V V 
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Til 


h \/+ \/v 



TABLE II: Diagram representation of the "bare" and the 
similarity-transformed second quantization operators. The 
black horizontal bars represent the cluster operators Ti and 
T<2, as displayed in Table U 

is not summed over. In practice, the circle prevents an 
accidental connection of the operator line, which would 
have introduced an erroneous Wick contraction when the 
spectroscopic factor diagrams are written down. 

The overall sign of a diagram is determined according 
to standard rules The negative sign of the second 

and third term in Eq. ([25)) is reflected in the at diagrams 
by the internal hole-lines that connect the small circle 
with the T operators. To determine the overall sign cor- 
rectly for the spectroscopic factor diagrams, a sequence 
of directed lines ending or starting in a small circular 
vertex must be counted as a loop. 

The diagrams shown in Tables Q] and [TT| are the basic 
building blocks for the computation of the spectroscopic 
factor. We compute the matrix elements of the overlap 
function as products of the components T M and ei- 
ther 4 or The only non-vanishing contributions to 
the spectroscopic factor come from the diagrams in which 
all directed lines can be connected. These diagrams and 
the corresponding algebraic interpretation are shown in 
Table IIIII We assume an implicit summation over re- 
peated indices. 

The computational cost of the spectroscopic factor di- 
agrams has the very gentle scaling, n^n^, so the cost is 
completely dominated by the CC and EOM calculations. 
In the case that \A) is the ground state of the closed-shell 
nucleus, we have rf = = r° h , and several diagrams van- 
ish. 



IV. RESULTS 

In this section we present our results for the calculation 
of the spectroscopic factor using ab initio coupled-cluster 
theory. We study the spectroscopic factor of nucleon re- 
moval from 16 by calculating the one-body overlap func- 
tions of 16 with the odd mass neighbors 15 and 15 N 
using the PR-EOM-CCSD approach to the ground and 
excited states of the A — 1 nuclei. The CCSD approxi- 
mation is used to calculate the ground state of 16 O. 

Our model space is spanned by oscillator states. We 
label the model space by the largest principal quantum 
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TABLE III: Diagram representation of the overlap expres- 
sions; repeated indices implies a summation. If the closed- 
shell system is in the ground state, diagrams involving either 
rf or rfj vanish. The individual components of these dia- 
grams are explained in Tables HI and HT1 



number N that is included in the single-particle (s.p.) 
basis, so that the maximum s.p. energy is En = (N + 
and the number of major oscillator shells is N + 
1. In Fig. Q] we show the convergence of the ground- 
state of 16 O with increasing size of the model space for 
a wide range of oscillator frequencies Hui, using Vi ow _fe 
with momentum cutoff A = 2.0fm~ 1 . In Fig. [5] we show 
the convergence of the ground-state energies of 15 O and 
15 N relative to the ground-state energy of 16 0. 

Our results shown in Figs. Q] and [2] show a weak de- 
pendence on hw in the largest model space. We esti- 
mate that our results for the ground state energies are 
converged within a few MeV in our largest model. We 
note that the CCSD ground-state for 16 is overbound 
by ~ 15MeV as compared to experiment. However, the 
energy difference between the ground-states of 15 and 
15 N is about ~ 4 MeV, which is very close to the ex- 
perimental value of 3.5 MeV. It thus seems that energy 
differences are better reproduced than absolute energies. 

The Huj dependence provides some information about 
how the finite size of the model space affects the solu- 
tions. For high values of hut, the model space includes 
high-momentum states beyond the momentum cutoff A 
of the interaction, but is not sufficiently extended in posi- 
tion space to accommodate a nucleus. For small values of 
Huj, the model space is sufficiently wide in position space 
for the extenson of the nucleus but does not contain suf- 
ficient high-momentum modes to resolve the cutoff A of 
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FIG. 1: (Color online) Convergence of the ground-state en- 
ergy (within CCSD) of ie O using a low-momentum poten- 
tial with cut-off A = 2.0fm~ for increasing model space size 
N — In + I and as a function of the oscillator spacing hw. 
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FIG. 2: (Color online) Convergence of the ground-state ener- 
gies (within PR-EOMCCSD) of 15 O and 15 N relative to the 
ground-state energy of 16 O using a low-momentum potential 
with cut-off A = 2.0fm~ 1 for increasing model space size and 
as a function of the oscillator spacing hw. 



the interaction. Close to the minimum, in the largest 
model spaces considered, a good compromise is realized. 

We also studied the energy levels using Vi ow -fe for var- 
ious momentum cut-offs in the range A = 1.6-2.2 fm _1 . 
The calculated ground-state energies for 16 0, 15 0, and 
15 N are sensitive to the cut-off, implying that induced 
three-body forces and short-ranged forces of higher rank 
would contribute significantly to the calculated energies. 

Let us turn to the spectroscopic factor for nucleon re- 
moval from 16 O. Figure [3] shows the spectroscopic fac- 
tor ([23)) 

SF(l/2-) = 5 1 1 |(/ = l,j = l/2) (30) 
for the removal of a proton with quantum numbers 



6 



i/2- 



from 16 O using a low- momentum interac- 
k with a cut-off A = 2.0fm _1 . Evidently, the 
spectroscopic factor is well converged and depends very 
weakly on the size of the model space and the oscillator 
frequency fno. It varies less than 1% over a wide range of 
oscillator frequencies. The spectroscopic factor SF(l/2~) 
for neutron removal from 16 O is almost identical to the 
SF(l/2~) for proton removal. Recall that isospin is ap- 
proximately conserved in light nuclei. 
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FIG. 3: (Color online) Spectroscopic factor SF(l/2~) for pro- 
ton removal from 16 O as a function of the oscillator spacing fko 
for different model spaces consisting of (iV + 1) oscillator shells 
and a low-momentum interaction with cutoff A = 2.0 fm - . 

The dependence on momentum cut-off A is displayed in 
Fig. 0J Note that the spectroscopic factor increases with 
decreasing cutoff. This is expected, since by lowering the 
cutoff the system becomes less correlated and the product 
state \4>o) becomes an increasingly good approximation, 
and the single-particle picture becomes more and more 
valid. Note also that isospin is approximately a good 
quantum number, as the spectroscopic factors for proton 
and neutron removal are almost identical. 

Let us also study the center-of-mass problem. The in- 
trinsic Hamiltonian ([TJ depends on the mass number A 
of the nucleus, and the calculation of the spectroscopic 
factor requires us to employ identical Hamiltonians for 
the nuclei with mass numbers A and A — 1. This con- 
stitutes dilemma, since no choice of actual value for the 
parameter A can satisfy the parent and daughter nuclei 
simultaneously. It is thus necessary to investigate how 
strongly the spectroscopic factor depends on this value. 
Figure [5] shows the spectroscopic factor (in a model space 
N = 4 for a momentum cutoff A = 2.0 fm -1 for different 
values of the mass number A of the intrinsic Hamiltonian. 
The dependence on A is very weak, and it is similar in 
size to the dependence on the parameters of the model 
space. 

For an intrinsic Hamiltonian, the coupled-cluster wave 
function of a closed-shell nucleus factorizes into an intrin- 
sic part and Gaussian for the center of mass of coordi- 
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FIG. 4: (Color online) Spectroscopic factor SF(l/2~) for neu- 
tron and proton removal as a function of the oscillator spacing 
hw for nucleon-nucleon interactions with different cutoffs in a 
model space with N = 6. 
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FIG. 5: (Color online) Spectroscopic factor SF(l/2~) for pro- 
ton removal from 16 O as a function of the oscillator spacing fku 
computed for different values of the mass number A employed 
in the intrinsic Hamiltonian {TJ. The model space consists of 
iV + 1 = 5 oscillator shells, and the momentum cutoff of the 
nucleon-nucleon interaction is A = 2.0 fm -1 . 



nate [13]. Following the procedure of Ref. [42j], we con- 
firmed that this factorization is present for the ground 
states of 15 O and 15 N in the largest model space we con- 
sidered. We found that this factorization even takes place 
if the value A = 16 for the mass number is employed 
in the intrinsic Hamiltonian (JTJ) for the computation of 
the nuclei 15 and 15 N. These results suggest that our 
approach to calculate spectroscopic factors within the 
coupled-cluster method is practically free of any center- 
of-mass contamination. 

So far, we focused on the spectroscopic factors for re- 
moval of a J 11 = l/2~ proton and neutron from 16 0. 
We finally also compute the spectroscopic factor for re- 
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moval of a J T = 3/2 ~ proton and a (deeply bound) 
= 1/2+ proton from 16 0. The result is shown in 
Fig. [6] for different model spaces. As before, the results 
are well converged with respect to the size of the model 
space, and display only a mild dependence on the os- 
cillator frequency. We find that the spectroscopic fac- 
tor SF(3/2") is similar in size to SF(l/2"). This is 
an interesting result. Barbieri and Dickhoff [25| also 
found in their computation of spectroscopic factors that 
SF(l/2~) « SF(3/2~) for nucleon removal from 16 0. As 
expected, the spectroscopic factor of the J 71 " = 1/2+ state 
is very small. The removal of a deeply bound = 1/2+ 
proton from 16 yields a highly excited state of 15 N that 
is a rather complex superposition of many n-particle- 
(ri + l)-hole states and thus has little overlap with a one- 
hole state. 
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FIG. 6: (Color online) Spectroscopic factors SF(3/2— ) and 
SF(l/2 + ) for proton removal from 16 as a function of the 
oscillator spacing hw. The employed model spaces have 7V+ 1 
oscillator shells, and the momentum cutoff of the nucleon- 
nucleon interaction is A = 2.0 fm~ . 



V. CONCLUSION AND OUTLOOK 

We have extended the coupled-cluster method for the 
computation of spectroscopic factors. To this purpose, 



we derived diagrammatic and algebraic expressions of 
the one-body overlap functions based on equation-of- 
motion methods for the ground and excited states of 
the closed-shell nucleus with mass number A and the 
neighboring nuclei with mass number A — 1. We imple- 
mented the equations in an uncoupled m-scheme, and 
presented proof-of-principle calculations of the spectro- 
scopic factor for proton and neutron removal from 16 0. 
The calculated spectroscopic factors are well converged 
in model spaces consisting of six oscillator shells for low- 
momentum nucleon-nucleon interactions. Within the 
coupled-cluster approach, the same intrinsic Hamiltonian 
has to be employed in the nuclei with mass numbers A 
and A — 1. We found that the spectroscopic factor is in- 
sensitive to the actual value of the mass number that is 
employed in the intrinsic Hamiltonian. 

We plan to implement the computation of the spec- 
troscopic factor also in a spherical formulation of nuclear 
coupled-cluster theory. This will allow us to employ much 
larger model spaces, and we plan to apply the techniques 
to physically interesting nuclei, such as 22 ' 24 0, 48,52 Ca, 
and 56 - 78 Ni. 
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